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Body movement monitoring device 

The present invention relates to a device for measuring body movement and more 
precisely for measuring gait parameters. 

5 The invention also relates to a method for using said device. 

Walking is one of the most common human physical activities. Evaluation of time and 
distance parameters during walking is helpful in assessing abnormal gait, to quantify 
improvement resulting from interventions, or to predict subsequent events such as falls. 
This evaluation is useful in several clinical situations such as: functional performance 
10 assessment after treatment or surgery such as hip and knee arthroplasty (Wall et al, 
1981), refining proper alignment and fit of external prosthesis or orthesis (Czerniecki and 
Gitter, 1996), fall risk assessment in elderly persons (Maki, 1997), or selecting the 
appropriate assistive device. 

A gait laboratory based on camera, walkway with implanted sensors (Hirokawa and 
15 Matsumara, 1987) or force-plates, and electromyography allows a complete gait analysis 
but require a dedicated laboratory. Although these techniques have been widely used for 
research purposes, their sophistication, time required for setting up the instrumentation 
and analyze the data, as well as cost, have hindered their use in clinical practice. In 
addition, these techniques require the subjects to walk in a pre-defined specific path, and 
20 assume that data measured from only a few steps are representative of usual gait 
performance. To avoid these limitations, ambulatory recording systems carried by the 
subject and allowing outdoor measurement have been developed, using new technology, 
such as a powerful microcontroler, miniature sensors, high capacity memory and small 
batteries. 

25 Ambulatory systems using foot switches or pressure sensors attached to the sole are used 
to monitor temporal parameters (Zhu et al., 1991; Abu-Faraj et al.,1997). These 
techniques provide satisfactory results for normal walking, but sensor attachment proves 
difficult when assessing subjects with abnormal gait. Even when correctly done, several 
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problems, such as shuffling gait, mechanical failure, or patient acceptance limit their 
applicability. Another system using insole with a matrix of multitude pressure sensors has 
a high reliability to evaluate the time of each phase of the gait cycle. Although this system 
provides additional information on foot pressure distribution during gait, it cannot be 
5 used for a long period of measurement, and sensors do not provide information during the 
swing phase. 

Recently, technical progress has made possible the realization of miniature kinematic 
sensors such as accelerometers and angular rate sensors with integrated conditioning and 
calibrating module (Sparks et al., 1998). In addition, due to their very low consumption, 

10 these sensors can be battery powered and are promising tools for long term ambulatory 
monitoring. The possibility of detecting simple parameters such as steps and cycle time 
from trunk or heel acceleration has already been demonstrated using miniature 
accelerometers (Evans et al, 1991; Aminian et al., 1995). More recently the efficiency of 
accelerometers to detect temporal features of walking has been shown (Aminian et al., 

15 1999). However, the acceleration signal is often biased by the gravitational acceleration 
and depends on the attachment site along the segment (e.g. limb). Fyfe et al. (2000) have 
proposed a motion analysis system, which measured the stride length and velocity of the 
foot. However, this system do not provided information related to the lower limbs 
rotation and temporal parameters. Since the actual nature of walking consists of lower 

20 limb rotation around joint articulation, the use of miniature angular rate sensors 
(gyroscopes) has proven to be an alternative technique for gait analysis (Tong and Granat, 
1999, Nene et al., 1999). However, no clear correspondence has been established between 
gait events and the patterns of angular velocity. This is probably because gait events (e.g. 
heel strike and toe-off ) are transitory signals that cannot be properly enhanced by 

25 traditional signal processing. 
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1 Summary of the invention 

One of the objectives of the present invention is to avoid the previous cited problems. 

It relates to a body movement monitoring unit as defined in claim 1 and to a 
corresponding method as defined in claim 7. 

5 In a preferred embodiment wavelet transform is used to overcome the previous discussed 
problems . 

What follows is a description of an ambulatory system according to the invention which 
is based on autonomous sensing units (ASU) for estimation of spatio-temporal parameters 
during long periods of walking. The method according to the invention is based on 

10 wavelet transform to compute the values of gait parameters from the angular velocity of 
lower limbs. The accuracy of measurements was assessed, using as a criterion standard 
the information provided by standard such as foot pressure sensors, forceplate and 
cameras. In order to estimate the accuracy of the method on a broad range of performance 
for each gait parameter, data were gathered from young (group A) elderly (groupe B) and 

15 subjects. The usefulness and the reliability of the proposed methods and device have been 
shown by performing gait analysis in three populations: patients with coxarthrosis (group 
C), patients with total hip arthroplasty (group D) and subjects age matched to the other 
two groups (group E). 

20 2 Methods 

2.1 Autonomous sensing units 

An autonomous sensing unit (ASU) was designed including a set of sensors, conditioning 
circuit (powering, amplifying and filtering), analog to digital converter, display, high 
capacity digital memory (64 Mbytes or more), rechargeable battery, and a microcontroller 
25 which controls the above components. The ASU was embedded in a small case and 
allowed saving signals detected by the sensors for a long period (fig 1). This 
configuration avoids the use of cable between sensors and a dattalogger and presents an 
easy way for attachment and comfort of the subject. The device allows also recording 
acceleration, angular rate in different axis of movement and also detecting temperature in 
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order to compensate kinematic signal (acceleration and angular rate) for temperature 
variation. Lower limbs movement during walking was measured using miniature angular 
rate sensors: gyroscopes. The gyroscope consists of a vibrating element coupled to a 
sensing element, acting as a Coriolis sensor. Coriolis effect is an apparent force that arises 
5 in a rotating reference frame and is proportional to the angular rate of rotation. Currently, 
Coriolis effect is used in several miniature gyroscopes such as piezoelectric gyroscope 
(Murata, ENC-03J) or MEMS gyroscope (Analog device, ADXRS300). Because of their 
low cunent consumption, these sensors can be battery powered and are therefore 
appropriate for ambulatory monitoring. 

10 Four ASUs including each at least one gyroscope were used. Using the conditioning 
electronic of the ASU, the signal from each gyroscope was amplified and low-pass 
filtered to remove the electronic noise of the gyroscope. A swiveling electromechanical 
system was designed to calibrate and characterize each sensor. Various angular velocities 
were applied to the gyroscope. An optical sensor measured the time corresponding to 

1 5 perform a rotation of 360° and allowed calibration of the gyroscope. The amplifier gain of 
each gyroscope was set to obtain a sensitivity of 4mV/deg/s. The ASUs were attached 
with a rubber band to each shank and on each thigh. Each ASU measured the angular rate 
rotation parallel to the mediolateral axis (perpendicular to the sagittal plane). The signals 
were digitized at a sampling rate of 200 Hz and stored on the memory of ASU. In order to 

20 synchronize the signal detected by each ASU and be able later to compare their values, 
before the measurement, all ASU were connected to the station unit which allowed 
starting all ASUs at the same time. During .the recording the ASUs were disconnected 
from the station and attached on lower limbs. At the end of the recording all ASUs were 
connected again to the station. The recording was stopped with the station unit. 

25 Afterward, the data saved on each ASU were transferred to a computer via the station unit 
and an USB port. The station unit has three major tasks: synchronizing all measurements, 
downloading recorded data to the computer and recharging the batteries of all ASU. 



30 



4 



WO 03/065891 



PCT/CH03/00091 



2.2 Validation with Footswitches 

Footswitches have been used as criterion standard to detect the gait events and validate 
the gyroscope measuring method The footswitches consisted of two pairs of force 
sensing resistors, FSR (Interlink, LU) placed under the right foot (fig.2): FSRhai beneath 
5 the heel (two sensors) and FSRt 0C beneath the big toe (two sensors). For each pair, the 
sensors were connected in series and the outputs were synchronized with those from the 
gyroscopes and recorded using a portable datalogger at a sampling rate of 200 Hz. During 
heelstrike and toe-off, FSRheci and FSR t0 c signals are subjected to abrupt changes. By 
thresholding the derivative of these two signals the exact time of heelstrike and toe-off 
10 was obtained. 

2.3 Validation with forceplate and camera based system 

Gait Analysis was performed by means of a motion analyzer (ELITE System, BTS SrL, 
Milano, Italy) equipped with 4 video-cameras recording at 100 Hz (Ferrigno G., Pedotti 

15 1985). One forceplate (Kistler 9661 A, 60 by 40 cm large) was embedded in the floor 
about one third down a 10 meters long walking pathway. Signals from the forceplate were 
collected at 500 Hz. The retroreflective markers were positioned on the lower limbs (fifth 
metatarsal head, lateral malleolus, lateral femoral epicondyle), on the pelvis (posterior 
iliac spines and lower edge of sacrum), and spine (posterior processes of seventh vertebra 

20 and point of maximum kyphosis). Joint kinematics was confined to only 
flexion/extension angles at the hip, knee, and ankle joints. 

The spatio-temporal parameters were computed on the basis of the heelstrike and toe-off 
(respectively initial and terminal contacts) of each foot on the ground. The forceplate was 
used to detect the initial contact of one foot. This was defined when the ground reaction 
25 force achieved a threshold of 5% of body weight. The next contact of the foot on the 
ground occurred in a non-instrumented area of the floor. Contact was defined when a 
similar markers' configuration was achieved (same distance from the floor, same 
displacement) as in the contact detected by the forceplate. The contralateral side was 
analyzed by looking at the distance of the foot from the ground. 

30 
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2 A Validation with treadmill 

In order to evaluate the accuracy of stride length and velocity, the following 
measurements were carried out. During treadmill walking, the value of treadmill speed 
has been considered as the actual mean walking velocity. A preliminary experiment in 
5 which various treadmill speeds were compared to the speed of an odometer placed on the 
treadmill belt confirmed that the treadmill speed is accurate. During ground (or "non- 
treadmill") walking, the mean velocity was evaluated by considering the time needed to 
cover 20m using a stopwatch. The actual mean stride length was calculated by 
multiplying the mean velocity by the gait cycle time obtained from footswitches. 

10 

2.5 Subjects 

Measurement using gyroscopes and FSR were taken from group A (9 young subjects, 
21± 2 yrs) and group B (1 1 elderly, 79 ± 8 yrs). Each young subject performed four trials 
including at least 20 gait cycles. The first three trials were performed on a treadmill at the 

15 preferred speed, under and over the preferred speed, respectively. The fourth trial was 
performed with a comfortable velocity over ground on a 30 m long walkway. The elderly 
performed only two trials on the walkway at their comfortable velocity. The study 
protocol was accepted by the institutional review board (Faculty of Medicine, University 
of Lausanne). Written informed consent was obtained from each subject. 

20 Measurements using gyroscopes, forceplate and Elite camera system were taken from 
group C (11 patients with unilateral coxarthrosis, age: 60±9 yrs), group D (8 patients 
with unilateral total hip prosthesis, age: 69±4 yrs) and group E (9 healthy subjects, 63±4 
yrs, weight: 63±9 kg). Each patient was asked to walk at a self-selected velocity down a 
walkway that included the forceplate. In order to reach good precision for the temporal 

25 parameters, each subject performed several trials (6 to 10 depending on patient condition, 
half with left foot touching the forceplate and half with right foot). At the end of the 
measurement, the exact time of toe-off and heelstrike extracted by both of the systems 
(gyroscopes and Elite/forceplate) were compared. 
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2.6 Wavelet analysis: Gait Temporal parameters estimation 

In order to compute the temporal parameters such as the duration of swing, single and 
double stances during a gait cycle, it is necessary and sufficient to determine for each leg 
the precise moments of heelstrike (when the foot first touches the floor) and toe-off (when 

5 it takes off) during that cycle. These events have distinctive signal features of shank 
angular velocity appearing as rather sharp negative peaks involving some medium and 
relative high frequencies (fig.3). Although the amplitude of these peaks vary according to 
various parameters such as subject's velocity or weight, or the presence of limping due to 
a painful articulation, they can always be localized as long as one knows where to look 

10 for it in time and frequency domains. To successfully achieve this task, we have 
decomposed the shank angular velocity into wavelet packages. Wavelet transformation is 
well adapted for gait events identification because it allows detection of a specified 
frequency at a specified time (e.g. toe-off and heel strike). In other words, toe-off and heel 
strikes events consist of combined features that can be well located in the time and 

15 frequency range. The methods of analysis need to be versatile enough to handle gait 
events that can be at opposite extremes in terms of their time-frequency location. 
Moreover, wavelets transformation allows the use of a suitable wavelet function (Coiflet 
wavelet), which is more similar to the pattern of heel strike and toe-off. 

A multi-resolution wavelet decomposition (Mallat, 1989) was used to enhance the toe-off 
20 and heel-strike event during walking. This technique consists in splitting a signal into 
high-scale (low-frequency components) called approximation and low-scale (high- 
frequency components) called detail. The decomposition is iterated, with successive 
approximations (or details) being decomposed in turn, so the original signal is broken 
down into many lower-resolution components. At each decomposition scale the number 
25 of samples in the time domain is decreased by throwing away every second sample (down 
sampling with a factor of '2'). By considering the original signal s(n) (shank angular 
velocity), the approximation of the signal at scale j=0 is A 2<) s which corresponds to the 

original signal s(n). At each new decomposition the approximation and detail decrease in 
resolution since at each scale 2 samples will be replaced by one. At scale j the A ,s 
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represents approximation of s(n) with a resolution of one sample for every 2 ] samples of 
the original signal. We have: 

s = A 2i s + D 2t s = A^s + D z2 s + D 2 ,s = ... = A 2 ,s + D 2j s + D 2 ^s... + D 2 ,s (1) 
Thus: 

5 A 2i s = A 2j . x s + D 2hl s (2) 

By using a suitable low-pass filter h 9 and a high-pass filter g, the approximate signal 
A lhx s and detail signal Z) 2/+1 .s can be further written as following: 

V>*= £*(2»-*Ki' (3) 

£g(2n-4V (4) 

10 Furthermore the signal can be reconstructed from approximate and detail signal using the 
following transform: 

A v s = 2 jT/fyi -2k\A lht s + 2 -2*)^.,* (5) 
Where, 

£(*) = /,(-*), (6) 

15 g(n) = g{-n) (7) 

The coefficients of h and g filter are associated with the shape of wavelet considered for 
the analysis. In this study, a decomposition into 10 scales with the Coiflet wavelet of 
order 5 has been used. 

The flowchart for the estimation of toe-off and heel-strikes are shown in fig.4. For clarity, 
20 let us consider only one shank. It is obvious that the algorithm will be applied in the same 

9 

way for both shanks. First an approximation of s(n) corresponding to the s a = ^D 2 /S was 
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obtained. With this approximation the original signal was only considered in the range of 
0.14 to 36Hz. Drift and high frequency movement artifacts were canceled in this way. 
Then two new approximations were obtained. The first one enhanced the heel strike 
component (A 2% s a ~A l9 s a ) and the second one enhanced the toe-off component 

5 (A 2 ,s a -A 29 s a ). For each of these approximations, the time corresponding to the global 

maximum values of the signal (ms) were detected. These values corresponded roughly to 
the time of midswing during a gait cycle (see fig.5), however their exact significance is 
not of interest to us. They will only be used as reference to select the intervals in which 
negative peaks reminiscent of heelstrikes and toe-offs are to be found. So we next looked 

10 for local minima inside intervals [ms + 0.25s, ms + 2s] and [ms -2s, ms -0.05s], 
assuming that a gait cycle is always less than 2s. The result was two different series of 
samples hs (heelstrikes) and to (toe-offs). The toe-off detection was refined by finding the 
local minima inside of the original signal s in the intervals [to, to + 0.075s]. Finally, at 
each cycle, the actual heelstrike and toe-off must be so that their difference is considered 

15 reasonable for walking. Thus, for each hs(k) only the nearest to(l) verifying the condition 
0.1s<hs(k)-to(l)<2.5s was saved. If no to(l) met this condition, h(k) was deleted and the 
next heelstrike was considered. Once we are in possession of hsi (left heel strike), to L (\eft 
toe-off), hs R (right heel strike) and to*(right toe-off), every temporal parameter of a gait 
cycle can be easily computed as percentage of gait cycle. These parameters are: 

20 - Duration of each gait cycle (measured at right or left leg) 

RGCT(k)= to R (k+l)~ to R (k) or LGCT(k)= to L (k+l)- to L (k) t l<k<N, (8) 

N: number of gait cycles 

- Left stance (time between left heel strike and left toe-off) 

LS(k)=to L (k)-hs L (k)) (9) 
25 - Right stance (time between right heel strike and right toe-off) 

RS(k)=to R (k)-hs R (k) (10) 

- Initial double support (time between right heel strike and left toe-ofif), known also as 
left double thrust support time (Wall et al., 1981)) 
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IDS(k)=tO L (k)~ hs R (k) (J J) 

- Terminal double support (time between left heel strike and right toe-off), known also as 
right double thrust support time (Wall et al.,1981) 

TDS(k)=to R (k+l)- hs L (k) (12) 

5 - Double support 

DS(k)=IDS(k)+TDS(k) (13) 

In order to reach a good precision for the temporal parameters, at least 20 gait cycles have 
been considered for analysis (N>20). For each trial, the value of temporal parameters was 
averaged over the N cycles. Statistical analysis is then performed to find the significance 
10 and accuracy of the parameters obtained by gyroscopes in comparison with the FSR based 
system. The root mean square error (RMSE) was used to analyze the average difference 
between estimated data by gyroscope and actual data obtained by FSR. 



^ {estimated - actual) 2 



RMSE = J n^bcrof samples 

" number of samples 



15 2.7 Temporal parameters estimation using IIR filtering 

The wavelet filter described above is a kind of bandpass filter. The low-pass part was 
used to "smooth" the signal a bit to increase the precision in locating the peaks and a 
high-pass part to get rid of drift. 

Another technique to cancel the drift would be the Infinite Impulse Response (IIR) filter. 
20 The IIR filter used has the following transfer function: 

= - IT 05) 

1-az 1 ' 

The value of a used during mentioned statistical study was 0.995. Using this filter the cut- 
off frequency of the filter will be 0.25Hz. To obtain precisely zero-phase distortion, the 
25 filter was applied on the input data twice. After filtering in the forward direction, the 
filtered sequence was reversed and run back through the filter. To detect heelstrike, the 



10 



WO 03/065891 



PCT/CH03/00091 



mentioned IIR filter was used and to detect toe-off, after applying the same IIR filter, a 
Finite Impulse Response (FIR) low-pass filter with pass-band attenuation of less than 
0.5dB and a cut-off frequency of 30Hz was used. To have zero-phase distortion the above 
forward and reverse filtering was used. This method of detection of toe-off and heelstrike 
5 was validated with 259 gait cycles recorded by both gyroscopes and force-plate. This 
filtering has the advantage to speed-up the calculation time. 

2.8 Gait model: spatial parameters estimation 

We propose a double segment gait model involving both shank and thigh. In this model, 
the swing phase is considered as a double pendulum model, while the stance phase is 

10 considered as an inverse double pendulum model. Consequently, knowledge of stance 
and swing phases obtained from the method described in 2.6 is necessary. Figure 6 (a) 
illustrates the rotation of the shanks and thighs during a complete gait cycle. The distance 
di+d2 corresponds to the part of stride length performing by the swing of right shank and 
right thigh. In addition, during this right swing phase, the body has moved forward by the 

15 distance d 3 involved by the rotation of left shank and left thigh throughout the 
contralateral stance phase. Stride length (SL) can then be measured as follows: 

SL=di+d 2 +d s (16) 

By considering a and p to be the angular rotation of the right thigh and right shank 
respectively, one can calculate the distance d=d|+d 2 from the trigonometric relations 
20 (fig.6.b). 

During the swing phase we have for each gait cycle k: 

d(k)=d,(k)+d 2 (k)= V(/ 2 +M l (*)) 2 + (/ 2 +M 2 (k)) 2 -(/ 2 + M l (*))(/ 2 +M 2 (*)).cosp(Ar) 

(17) 

where: 
25 //: thigh length 
l 2 : shank length 



sin P(fc) (18) 
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sin (3 (k) ( ]9 j 
wi th: rW = ^^ and and -^2(1- ««<*)) 

The angles a and P were estimated by integration of the angular rate rotations of the right 
thigh (t R ) and right shank (s R ) as follow: 

<*(*) = \t R {t)dt P(*)= J^W-A (27; 

'M*> w«(*) 

During the stance phase we have: 

d(k)=d 3 (k)=J{l x + M 2 (£)) 2 + (/, + A/, (£)) 2 - (/, + Af 2 (£)) * (/, + M 1 (*)).cosa(Ar) f»; 
where: 

sinP(*) ^ 

w= sinK*) d , ()t) 

smp(*) ^ 

with:y(*)= g " PW and <p{k) = ff ~ 2a W + P<*> 3 nd rf'(Jt) = / 2 (*)V2(l-cosp(*)) 
and: 

<*(*) = Jf*(0'* P(*) = (26) 

For each gait cycle k, the stride length (SL) and stride velocity (SV) were obtained as 
follow: 

SL(k)=d,(k)+d 2 (k)+d } (k) (27) 

SV(k)= SL( - k) (28) 
l/ RGCT(k) 1 ' 



12 



WO 03/065891 



PCT/CH03/00091 



where RGCT(k) was obtained from (8). 

If necessary, we can use only three gyroscopes instead of four, by assuming the same 
length for left and right steps, one can consider that the rotation during left stance phase is 
equal to rotation during right stance phase. In this case (26) change to: 

5 a(*)= \t R {t)dt P(*)= \s R {t)dt (29) 

hs R {k) hs R (k) 



2,9 Comparison with references systems 

Toe-off and heelstrike events were estimated using FSR (footswitches) in group A and B 
and using forceplate and Elite camera based systems in groups C, D and E. At the end of 
10 the measurement, the exact time of these events were compared to those estimated based 
gyroscopes. 

In order to evaluate the accuracy of stride length and velocity given by (27) and (28), the 
following measurements were carried out. For group A treadmill will used, the value of 
treadmill speed has been considered as the actual mean walking velocity (V a ). For group 
15 B only overground walking was performed. The mean velocity was evaluated by 
considering the time needed to cover 20m using a stopwatch. The actual mean stride 
length (L a ) was calculated as: 

L 0 =V a GCT a (30) 

where GCT a is the mean gait cycle time obtained from footswitches. The accuracy of 
20 stride length and velocity was evaluated by comparing L a and V a with the mean of SL and 
SV{SL m and SV m ) over 20m. 

In group C, D and E stride-length and velocity (SL a and SV a ) were obtained for each cycle 
based on Elite and forceplate systems as explained earlier and compare to our estimated 
values (SL m and SV m ). 

25 
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3 Results 

3.1 Temporal parameters 

The correspondence between temporal events detected by FSR and the gyroscope pattern 
is shown in fig. 5 for five consecutive gait cycles. Using the algorithm of fig. 4, left and 
5 right heel strikes (respectively toe-off) were detected from shank gyroscopes (fig. 7). Then 
based on relation (8) to (13), the mean value (over the N cycles) of left and right gait 
cycle time (LGCT^ RGCT m \ left and right stance (LS m and RS m ), initial and terminal 
double stance (IDS m , TDS m DS m ) were computed for each trial on treadmill and ground. 
Simultaneously, the periods of the actual gait cycle (GCT a ) and right stance (RS a ) were 
1 0 detected from FSR hee i and FSR t0 e. 

Figure 8 compares stance period (RS m and RS a ) and gait cycle times values (RGCT m and 
GCT a ) obtained from the gyroscopes to those of FSR from a trial performed by the group 
A (young subjects: treadmill and ground) and group B (elderly during overground 
walking). A good agreement (r>0.99) was found. The difference between RS m and RS a 

15 was not significant (N=58, p>0.80) with a mean standard error RMSE = 23ms. No 
significant change was observed between RGCT m and GCT Q (N=58, p>0.97). The RMSE 
between RGCT m and GCT Q was 8ms (0.5%). The 95% confidence interval for the 
difference between heelstrike detection by FSR and gyroscope was [7ms, 13ms]. 
Heelstrike detection based on the gyroscopes occured later (10ms in average) than FSR 

20 (underestimation). The 95% confidence interval for the difference between toe-off 
detection by FSR and the gyroscopes was [-5ms, 4ms]. Therefore there was no significant 
difference in toe-off detection. The RMSE between LGCT m and RGCT m was 15ms 
(r>0.99, 1 .3%), which shows the validity of the method in providing the same time for 
both left and right gait cycles. 

25 The mean values of right and left stance (RS m and LS m ) as well as double supports (TDS m , 
IDSm, DSm) estimated by the gyroscopes are shown in table 1. In general, higher values 
are observed for elderly subjects (group B). 
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Subjects 


RS m % 


LSm.% 


IDS m % 


TDS m ,% 


DS m ,% 


Group A(Young) 


61 ±2 


62 ±2 


12 ±2 


12 ±2 


23 ±3 


Group B(Elderly) 


64 ±3 


65 ±2 


14 ±3 


15 ±2 


29 ±4 



Table 1- The averaged values of right stance {RS m ), left stance {LS m \ initial double stance 
5 (IDS m \ terminal double stance (TDS m ) and double stance (DS m ) in the percentage of gait 
cycle obtained for young (group A) and elderly subjects (group B). 

A total of 106, 80 and 83 gait cycles were obtained from groups C, D and E, respectively. 
A good correspondence was observed between gyroscopes and the Elite system. By 

10 considering all the gait cycles (N=269), a systematic difference of -14ms was observed 
between the actual initial contact (from the forceplate) and the corresponding value 
detected by gyroscopes. This difference reached 4ms for terminal contact. Gyroscopes 
detect slightly later the heelstrike of the foot, while toe-off detection based on gyroscopes 
occurred a little earlier than with the forceplate. The same trends were observed within 

15 each group C, D and E. Assuming this difference as a systematic error, the values of toe- 
off and heelstrike were corrected by considering these differences. Students t-tests applied 
on temporal parameters showed no significant difference (p>0.99) between the actual 
values of gait cycle time and stance period, as compared to those estimated by 
gyroscopes. Table 2 illustrates for group C, D and E, the RMSE and the overall 

20 coefficient of correlation versus the actual values. The error ranged between 2.7% and 
4.6% (3.7% when considering all subjects). 

3.2 Spatial parameters 

In group A and B, for each trial the mean value of stride length and velocity (SL m and 
25 SV m ) were estimated by finding the average of stride length (SL) and stride velocity SV 
(relation 27 and 28) over the 20m of walking. Fig. 9 compares actual L Q and V Q to SL m and 
SV m . No significant change was found between L a and SL m {N=56, p>0.52) and neither 
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between V a and SV m (p>0.66). The coefficient of correlation is high in both cases 
(rXX96). Actual velocity and stride length can be estimated from the measured values as 
follow: 

V a =L02 SV m -0.04 (31) 
5 L a =0.97SL m (32 ) 



The RMSE for velocity estimation is 0.06m/s (6.7%) and that of stride length estimation 
is 0.07m (7.2%). 

Stride length and velocity had lower averaged values for the elderly (group B) compared 
1 0 to young subjects (group A). 

In groups C, D and E, stride length and velocity, shows also a linear relationship between 
actual (SV at SL a ) and gyroscopes values (SV mt SLJ. When including all subjects in these 
3 groups, we obtained the following linear regression: 

15 SV a =0.931 SV m +0.006 
SL a =J.172SL m -0.214 



(33) 
(34) 



When considering each group C, D, E separately, the slope and intercept were also close 
to the value of equation 31 and 32. Based on these linear approximations, the values of 
20 SV m and SL m were modified to be closer to the actual values (replacing the measured 
value by its linear approximation). The coefficient of correlation (r) and RMSE for spatial 
parameters are shown in table 2. The error ranged between 4.8% and 7.5% (overall error 
of 7%). No significant difference were observed between the actual values and those 
obtained from gyroscopes(p>0.99). 

25 Figure 10 compares the ranges of thigh and shank rotation estimated from gyroscopes 
with the corresponding values of the Elite system. A very good agreement (r=0.97) 
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between the two systems can be observed. The overall RMSE was less than 5 degrees 
(Table 2). 





RMSE 


r 


Group 


C 


D 


E 


Overall 




Initial contact 


16ms 


9ms 


11ms 


13ms 


>0.99 


Terminal contact 


23ms 


34ms 


24ms 


28ms 


>0.99 


Gait cycle 


38ms 


52ms 


30ms 


38ms 


0.93 


Stance 


25ms 


33ms 


28ms 


28ms 


0.92 


Stride velocity 


0.07m/s 


0.06m/s 


0.08m/s 


0.07m/s 


0.88 


Stride length 


0.08m 


0.05m 


0.08m 


0.08m 


0.82 


Thigh rotation range 


4deg 


7deg 


5deg 


5deg 


0.82 


Shank rotation range 


3deg 


2deg 


2deg 


3deg 


0.94 



5 Table 2. Error (RMSE) for gait events (toe-off and heelstrike, gait cycle, stance time) and 
spatial parameters (stride length and velocity, range of thigh and shank rotation) obtained 
for groups C, D, E and by considering all subjects of these three groups (overall). The 
overall correlations (r) show the good agreement between our proposed method and 
reference system. 

10 

Figure 1 1 shows typical results obtained during 2 minutes of treadmill walking where the 
velocity was fixed at 1.11 m/s (4 km/h). The subject was asked to walk at his normal 
pace, to decrease his pace then increase his pace. All parameters were obtained at each 
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gait cycle from gyroscopes. As it is illustrated, the estimated velocity (SV) is almost 
constant {mean-\.\A m/s) with a coefficient of variation of 4%. The change of gait cycle 
(RGCT) provides a corresponding change of stride length (SZ) in order to maintain the 
velocity constant. The corresponding coefficient of variation of SL and RGCT is 15%. 
5 During the same period there is no significant variation of left and right stance 
(coefficient of variation less than 3%) while the double support time (DS) has a 
coefficient of variation of 16%. 

Toe-off and heelstrike detection based on IIR and FIR filtering are compared favorably 
with wavelet techniques. In 259 gait-cycles the measured signals of gyroscopes were 
10 analyzed twice, once using the wavelet based filter and once using an IIR + FIR filter 
method. Table 3 compares the obtained results by using forceplate as reference. 





Mean (ms) 


Std (ms) 


Wavelet 


Toe-off 


4.6 


29.2 


Heel-strike 


-13.3 


14.0 


IIR+FIR 


Toe-off 


-2.9 


28.0 


Heel-strike 


-13.6 


13.2 



Table 3. Toe-off and heelstrike detection by using the two proposed method of filtering: 
1 5 Wavelet filtering and IIR+FIR. Mean and standard deviation are expressed for 259 gait 
cycles. 

4 Discussion 

We have proposed a new algorithm based on wavelet transformation to detect toe-off and 
20 heelstrike from shank angular velocity. A high correlation with actual gait events detected 
by standard reference systems (foot pressure sensors, forceplate and camera) was found. 
In addition, no significant error was observed for temporal events detection. The slight 
delays that can be observed between reference systems and gyroscopes are in the range of 
resolution of the used systems. In addition, the obtained values correspond to the 
25 established values for young and elderly populations which is around 60% for the 
duration of stance and 20% for the duration of double support ( Rose et al., 1994, Maki, 
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1 997). Few studies have been reported about temporal gait parameters 1 detection based on 
gyroscope. Pappas et al. (1999) used a gyroscope on the heel in conjunction with FSR to 
detect the precise time of heel-off. Tong and Granat (1999) used a gyroscope on the shank 
and thigh. They found that the pattern on the shank has two minima, one occurs at foot 
5 flat and the other occurs at toe off. This last result is consistent with ours. In contrast to 
their results, we found that the first peak after wavelet enhancement corresponds to heel 
strike and not foot flat. It must be noted that only 2 subjects were involved in their study 
and the sampling rate was 50Hz. Toe-off and heelstrike detection were based on the 
choice of a threshold on the FSR signal. It is clear that the value of the threshold can 
10 modify the results. Moreover, only a few steps were considered for analysis which can not 
be considered as statistically significant. 

The method presented in this study compares favorably with other methods used to 
identify gait events that require more complex instrumentation and a laboratory setting. 
Hreljac and Marshall (2000) used four motion analysis cameras (60Hz), involving 12 

15 trials of two healthy subjects and found an absolute average error of 4.7ms and 5.6ms for 
heelstrike and toe-off, respectively, considering information from a forceplate as the "gold 
standard". Stanhope et al. (1990) using a 50Hz camera collection system, reported an 
error greater than 20 ms in over 20% of the cases. Using accelerometers on the thigh and 
considering FSR as reference, Aminian et al.(1999) found the 95% confidence interval of 

20 stance period equal to [-20ms, 20ms]. In the present study, which involve young and 
elderly subjects the error of toe-off and heelstrike detection is less that 5ms (95%), with a 
systematic delay for heelstrike time of 10ms. 

Gait events can also be explained by looking at the particular pattern of the shank and 
thigh angular velocities (fig. 5). Shank angular velocity (signal s R ) is a bipolar signal. 

25 Considering as positive the counter-clockwise rotation, one can observe that during the 
swing period s R is rather positive while in stance phase s R is almost negative. Just prior to 
heelstrike, s R becomes negative. At heelstrike, there is a flexion (Vaughan et al., 1992) of 
the knee shown by a negative s R and positive t R . The time of heelstrike is also 
characterized by a change of slope sign of s R . This is expected since before heelstrike, the 

30 shank decelerates (negative slope) to stop and after heelstrike there is an acceleration 
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(positive slope) due to rotation around the ankle. After heelstrike, there is a clockwise 
rotation of the shank around the ankle involving a negative sr. In mid-stance, the knee 
becomes straight leading to a slight angular velocity, which could become positive. Then 
the shank continues its clockwise rotation around the ankle and accelerates its rotation at 
5 the toe-off (positive slope of sr). The rotation is still clockwise (sr negative) at the 
beginning of the swing phase until the knee reaches its maximum flexion (t R is positive 
and sr negative), then the knee comes into extension (s R positive and t R positive) moving 
the shank from backward to forward. In this way, shank and thigh angular velocities also 
give remarkable insight into the recognition of the periods of knee flexion and extension 
10 during gait cycle. 

In order to measure spatial parameters such as stride length, we have used the change of 
lower limb angles. Miyazaki (1997) has shown that for the estimation of the stride length 
from the integral of angular velocity, the relative change of angular rotation is sufficient 
and therefore there is no need of initial angular position. However this method consider 

15 thigh and shank as a unique segment and the error reached 15%. We propose a double 
segment gait model involving both shank and thigh, which consider a double pendulum 
model during the swing phase and the inverse double pendulum model during the stance 
phase. The identification of swing and stance phases was therefore necessary, and was 
determined by toe-off and heelstrike estimation. The estimation error for velocity and 

20 stride length is around 7%. This error is principally due to the inaccuracy of the 
gyroscopes which induce some error on a and p in (20) and (25). This error could 
increase if the period of integration increases considerably (very slow speed). In this case 
a new calibration of (29) and (30) is necessary. Other factors such as the difficulty to 
measure the exact shank and thigh length, the misalignment of the gyroscope axis 

25 regarding the mediolateral axis and the error due to gait phase detection are also to be 
considered. The low values of the intercept (close to 0) and the slope close to 1 in (29) 
and (30) confirm the appropriateness of the proposed gait model. 

There are many advantages to use gyroscopes instead of other kinematic sensors such as 
accelerometers. First, unlike the accelerometer, the gyroscope can be attached to 
30 anywhere of any body segment as long as its axis is parallel to the mediolateral axis: the 
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angular rotation is still the same along this segment. Tong and Granat (1999) have shown 
that the signals from different gyroscopes at different attachment site are almost identical. 
Second, the angular rate signal is less noisy than acceleration since acceleration is the 
derivative of velocity and involves higher frequency components. Third, accelerometry 
5 provides essentially temporal gait parameters, using singular peak enhancement. In 
contrast, rotation angles can be estimated from angular velocity by simple integration, 
and provides spatial gait parameters in addition to temporal parameters. Recently, the 
possibility to estimate knee angle (Williamson et al, 2001) and trunk tilt (Najafi et al, 
200 l a ) has also been shown. Finally, there is no influence of gravity acceleration on the 

10 measured signal. Gyroscopes however, do have some weaknesses. First, gyroscope is 
more delicate to use than accelerometer and can be more sensitive to temperature. In 
addition, powerful signal processing and filtering are necessary to cancel drift and 
artefact in the signal. Nevertheless, this filtering was done in this study by using wavelet 
transformation. The use of wavelet was particularly important to ensure a good time 

1 5 resolution in finding gait events such as toe-off and heelstrike. 

Overall, the autonomous measuring devices (ASU) are light (each 80 gr), portable, 
wireless, low cost and involves no discomfort for the patient who carries it for a long 
period of time. The ASU is battery powered with a lifetime of 20 hours. It can collect 
data up to 64 Mbytes corresponding to a period of 42 hours (at 200Hz sampling rate). 
20 The number of sensor to be used depends on application. Two gyroscopes can be used on 
one side (e.g. right shank and right thigh) where the estimation all temporal parameters is 
not desired. By using gyroscopes on both sides the accuracy of the system will increase 
and higher number of gait parameters will be estimated. 

Finally, the feasibility of the measurement in pathologic cases, where extremely low 
25 movements or, brisk contractions of muscles can affect the quality of the signals has been 
shown. One major feature of the proposed method is that we have nearly the same 
accuracy for all three groups (coxarthrosis, prosthetic hip and control). 



21 



WO 03/065891 



PCT/CH03/00091 



5 Conclusion 

The proposed method appears a promising monitoring tool for several purposes. First, it 
allows measurements of gait features during a long period of walking and thus supplies 
the stride-to-stride variability of gait as illustrated in fig. 10. In addition, because of its 
5 portability, this system can be used in other settings than a gait laboratory and therefore 
provides information that is more likely to reflect the actual performance of the subjects. 

It can be used in many clinical applications such as outcome evaluation after knee and hip 
replacement, or external prosthesis adjustment for amputees. In elderly subjects, this 
system can also be proposed as a diagnostic tool for abnormal gait analysis, as a predictor 
10 tool for fall risk estimation (Najafi et al., 2001 b ), or as a monitoring tool to assess 
progress through rehabilitation. During this study the system proved to have a high 
acceptability by elderly subjects. 



15 Brief description of the figures 

Fig.l Autonomous Sensing Unit (ASU) with their electronic and sensors. ASU includes 
one or several gyroscopes allowing the angular rate measurement of body segment. 
Before starting measurement all ASUs are linked to a station for synchronization. 
Afterward, they are disconnected from the station and record segment's movement. At the 
20 end of measurement all ASU are again connected to the station. The recorded data are 
transferred to the PC via station. The station allows also recharging the battery of each 
ASU. 

Fig.2 Footswitches (FSR) attachment beneath the heel (two sensors) and beneath the big 
toe (two sensors). 

25 Fig.3 Shank angular velocity showing the presence of peaks during the toe-off and 
heelstrike. 
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Fig.4 Algorithm flowchart for temporal parameters estimation from shank angular rate 
signal (s). 

Fig.5 FSR signal from heel (FSRheei) and toe (FSR toe ) 9 angular rate signal at right shank 
(s R ) and right thigh (t R ). Heelstrike (0) and toe-off (o) detected by FSR and s R are shown 
5 on each signal during 5 typical gait cycles. m s correspond approximately to the time of 
mid-swing. 

Fig. 6 (a) Body segments model during a gait cycle started by right toe-off Qor) and ended 
by right heelstrike (As*), (b) Shank and thigh configuration during stance and swing 
phase. Thigh and shank length (/; and l 2 respectively) are measured on each subject, a and 
10 p obtained from angular rate signals (t H and s R ) allow to estimate the distance d during 
stance and swing phases. 

Fig.7 Temporal parameters estimation from heelstrikes (o) and toe-offs (□) of each leg 
based on relation 8 to 13 in the text. 

Fig. 8 (a) Right stance duration, in s, obtained from gyroscope (RS m ) compared with that 
15 obtained from pressure sensors (RS a ), r>0.99. (b) Gait cycle time from gyroscope 

(RGCT m ) compared to actual gait cycle time obtained from FRS (GCT a ), r>0.99. (•) 
young, (o) elderly subjects. 

Fig.9 Correspondence between actual stride length (L a ) and velocity (V a )and estimated 
values from sKank and thigh gyroscopes (SL„u SV m ) for young (•) and elderly subjects 

20 Fig. 10. Range of thigh and shank rotation obtained from the proposed method 
(Measured) and Elite (Actual). A: group C, O: group D, □: Group E 

Fig. 11 Stride velocity (SV), stride length (££,), gait cycle time(/?GC7), right and left 
stance (RS and LS) and double support (DS) in the percent of gait cycle obtained from a 
young subject during 2 minutes of treadmill walking using the proposed method. The 
25 velocity of treadmill was kept at l.llm/s while the subject changed his cadence. All 
parameters are estimated at each gait cycle allowing to analyse the stride to stride 
variability. 
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Claims 

1 . An autonomous measuring unit comprising : 

a. a set of sensors, 

b. conditioning means for deriving information from said sensor, 

c. display means for displaying said information to an operator, 
characterized by the fact that it furthermore comprises means for recording 
kinematic parameters of body segment. 

2. Autonomous measuring unit according to claim 1 wherein at least one sensor is an 
angular rate sensor designed for measuring rotation of the segment on which it is 
attached. 

3. Autonomous measuring unit according to any of the previous claims furthermore 
comprising means for measuring and recording temperature. 

4. Body movement monitoring system comprising a set of autonomous measuring 
units as defined in any of the previous claims, wherein the said autonomous 
measuring units are designed to be attached on lower limbs. 

5. Body movement monitoring system according to claim 4 furthermore comprising 
means for determining the time of initial and terminal contact (heelstrike and toe- 
off) of each foot, the range of rotation of each shank, each thigh and each knee, as 
well as the stride length and stride velocity during walking. 

6. Body monitoring system according to claim 5 wherein the means as defined in 
claim 5 include processing methods to detect the local minimum negative peaks of 
the shank angular velocity in order to find initial and terminal contact of the foot, 
and a double segment gait model involving both shank and thigh in order to 
estimate angular rotation of lower limbs, stride length and stride velocity. 
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7. Method for the use of a system according to anyone of claims 4 to 6 wherein all 
autonomous measuring units are synchronized by connecting to a station before 
the beginning of measurement and at the end of measurement. 

8. Method according to claim 7 wherein the data saved on each unit are transferred 
to the computer, and wherein the station maintains the charge of each unit and 
tests the state of the sensors and the state of any other essential elements such as 
the measuring circuits, the state of the memory, the battery and the display of each 
unit. 

9. Method according to claim 7 or 8 wherein the spatio-temporal parameters of gait 
are estimated during each gait cycle (two steps), providing the estimation of mean 
value and inter-cycle variation of these parameters during a walking trial. 
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